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SUMMARY 


The geometry of general three-dimensional bodies is generated from coordi- 
nates of points in several cross sections. Since these points may not be 
smooth, they are divided into segments and general conic sections are curve fit 
in a least-squares sense to each segment of a cross section. The conic section* 
are then blended in the longitudinal direction through longitudinal curves 
which may be used to define the conic sections in the cross-sectional pleines. 
Both the cross-sectional and longitudinal curves may be modified by specifying 
particular segments as straight lines or specifying slopes at selected points. 
Slopes may be continuous or discontinuous and finite or infinite. 

This method was used to siirface fit a 70° slab delta wing and the HL-10 
Lifting Body. The results for the delta wing were very close to the exact ge- 
ometry. Although there is no exact solution for the lifting body, the surface 
fit generated a smooth surface with cross-sectional planes very close to pre- 
scribed coordinate points. 


INTRODUCTION 


Many disciplines require a mathematical description of three-dimensional 
surfaces which cannot be represented by simple mathematical expressions. The 
location and slope of points on a body are needed in the analysis of structures 
and inviscid flow fields. Viscous flow-field analyses require the body curva- 
ture in addition to location and slopes. The geometrical properties of simple 
shapes like spheres , cones , ellipsoids , and paraboloids can be described by re- 
latively simple mathematical equations. However, many configurations of inter- 
est today, such as the space shuttle, are complex three-dimensional shapes 
whose geometry cannot be described very easily. In many cases all the informa- 
tion that is known is a drawing with a plan view, a side view, and several 
cross sections of the body. Sometimes models of a vehicle are available and 
coordinate positions can be accurately measured on them. On the other hand, 
slopes cannot be accurately measured on a drawing or model, and the measured 
radii of curvature are even less accurate. This report develops a method for 
surface fitting mathematical relations to complex three-dimensional bodies. 

The method yields accirrate coordinates and slopes and reasonably accurate radii 
of curvature at any position on the body. 

Previous approaches to surface fitting three-dimensional bodies generally 
divided the surface into "patches" and represented each patch by flat surfaces, 
cubic or higher order polynomials, or conic sections (e.g. see refs. 1-5). 
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Each of these methods has iindesirable features in the forms used previously. 
Flat surfaces are completely unacceptable if radii of curvature or continuous 
slopes are needed. Cubic or higher order polynomials often lead to unwanted 
wiggles and bulges since they allow points of inflection. In ref. 1, 

DeJarnette applied the method of double splines, which used bi-cubic interpola- 
tion, to surface fit the coordinates of points in several cross-sectional 
planes of three-dimensional bodies. However, it was found that bulges and/or 
dimples occiirred in these surfaces, particularly when the thickness was much 
smaller than the span. 

Reference 2 approximated the shape of a space shuttle with elliptical 
cross sections with different ellipticity on the windward and leeward sides. 
Cubic polynomials were used to define segments of the plan and thickness dis- 
tributions. The coefficients of these polynomials were chosen to make the 
slopes continuous across bo\andaries of the longitudinal segments. However, 
points of inflections were found to occ\ir inside the segments and thus gave un- 
desirable bvunps and wiggles, which in turn significantly affected the surface 
pressure and heating distributions calculated with this geometry. 

For sometime conic sections have been used to describe segments of air- 
craft contours (ref. 3). Both longitudinal and transverse contours were repre- 
sented by conic sections, but the slopes at the ends of each segment had to be 
measured. As mentioned earlier, slopes are difficult to measiire accurately. 

On the other hand, conic sections cannot have inflection points, and this fea- 
ture is useful in eliminating \uidesirable wiggles. 

Coons (ref. U) developed a technique to describe three-dimensional sur- 
faces by using blending functions to blend the surface between the boundary 
ciirves of each patch. A major difficulty of applying Coons' method is that the 
user must supply the coordinates, slopes, and twists (cross derivatives) at all 
fo\jr corners of each patch. This information is generally difficult to deter- 
mine, particularly the cross derivatives (See ref. 5). Craidon (ref. 6) used 
this method with the twists arbitrarily chosen to be zero to surface fit three- 
dimensional geometries. 

The method presented here uses data points in cross-sectional planes to 
ciirve fit segments of general conic sections in a least-squares sense. The 
conic sections are then blended in the longitudinal direction by use of coordi- 
nate points which define the conic sections in the cross-sectional planes. 

This technique has the advantage of allowing the user to continually modify the 
cross-sectional curves and the longitudinal curves until the body shape has the 
desired features. Also, discontinuous slopes in both the circumferential and 
longitudinal directions may be specified. 


SYMBOLS 


a b . parameters defined by Eqs. (B5) - (B16) 
ni , i » j 


A^ thru A^ 


coefficients of conic section given by Eq. (l8) in global coordinates 
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d 

P 

G 


pq 


K 


J 





coefficients of conic section given by Eq. (l) in local coordinates 

constant vector in Eq. (l9) 
matrix used in Eq. (19) 

data point number of first control point in segment J 
slopes defined by Eqs. (5) and (6) 


N number of segments in a cross section 

P , Q terms defined by Eqs. (A2) and (A3) 

J J 

r ,, residual in Eqs. (Bl) and (B2) 

J ^ 

R term defined by Eq. (A5) 

J 

X, y, z Cartesian coordinates, see figure 2 

y, z local coordinates, see figure 3 

“.J S., Yj coefficients defined by Eqs. (15), (l6), and (17) 

J J J 


A9 


J 


slope of segment j , see figure 5 




Subscripts : 


h intermediate point in a cross-sectional segment 

J segment number in a cross section 

k data point number in a cross section 


ANALYSIS 


In a typical application, the geometry of a body must be determined from a 
model or a three-view drawing. Generally several cross sections are obtainable 
from the model or drawing, and the coordinates of data points on the boundary 
of these cross sections can be measured or calculated. A three-dimensional 
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surface must then be fit through all the cross sections. In the method pre- 
sented here, the data points in each cross section are divided into segments 
and general conic sections are fit to the data points in each segment as shown 
in figiire 1. A three-dimensional surface is then generated by "blending" the 
cross-sectional ciirves in the longitudinal direction. Consider first the curve- 
fitting of data points in a cross-sectional plane. 


Curve-Fit in a Cross-Sectional Plane 

The data points in a cross-sectional plane are generally not completely 
smooth, and in these cases a smooth curve cannot be made to pass through every 
data point. Therefore, the data points are divided into segments, and a gener- 
al conic section is fit in a least-squares sense through the data points in 
that segment. The conic section is constrained to go through the control 
points (end points of a segment) and aiso have a continuous slope at each con- 
trol point unless otherwise specified (See figure l). 

A three-dimensional coordinate system x, y, z is used with x in the longi- 
tudinal direction and x = constant is a cross-sectional plane (See figure 2). 

Let J denote the segment n\jmber in a cross-sectional plane as shown in figure 1, 
with the first segment starting on the positive y - axis and J increasing clock- 
wise. It is convenient to use a local coordinate system y,z for each segment 
with the origin on the first control point and the positive y - axis passing 
through the control point at the other end of the segment (See figure 3). In 
this local coordinate system it is easy to investigate possibilities of complex 
roots and interpret the coefficients of a general conic section geometrically, 
whereas it is difficult to interpret them geometrically in the global coordi- 
nates y,z. It is assumed that the local coordinate z is single-valued in each 
segment . 

The equation for a general conic section for the j segment is given by 
ref. 7 as 


Ajy^ + B^yz + C^z^ + D^y + E^z + Fj = 0 . (l) 

Only five of the six coefficients are independent since the equation may be 
divided by any non-zero coefficient. The constraints that the conic sections 
pass through the two control points 

y = 0, z = 0 and y = y.» z = 0 yields 

J 

Fj = 0 (2) 

< 3 ) 

By differentiating Eq. (l) with respect to y, the slope in the local coordinates 
is given by 
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( 4 ) 


dz 

dy 




The slopes at the control points are generally not known, so define 


dz 

“"j ” i<iyJ 


at y = 0, z = 0 


and 




dz 


l<iyj 


at y = yj, z = 0 


By using these two equations in Eq. {k) it follows that 
and 




_1 _ 1 _ 
1“J " "JJ 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


These equations show that the end slopes do not affect the coefficient Cj. On 
the other hand, the product AjCj determines the nature of the general conic sec- 
tion (ref. j). If the discriminant (B? - Ua.C.) = 0 the conic is a parabola, 
if the discriminant < 0 the conic is sin ellipse, and if the discriminate > 0 
the conic is hyperbola. Of particular interest is the possibility of complex 
roots when solving for z as a function of y in the region of interest. It is 
shown in Appendix A that z will not be complex in the region 0 ^ y < y if 


A C > 

J J - 






( 9 ) 


For prescribed slopes mj , n. figure U shows how the product A.Cj affects a 
conic section. 


Unless a slope is specified at a control point, the method used here con- 
strains the slope in global coordinates (y,z) to be continuous at a control 
point. This makes the conic section in one segment dependent on data points 
in other segments as well as its own. From figure 5 it can be seen that con- 
tinuity of slope at control point j requires 


tan ^ ^^41 ~ ^ • 

J J J 

This equation can be expanded and rearranged into the form 

A. 


. 1-1 
TaTS" 


J-1 


(A/m) sin A0 + A cos A6. 

J J J- 

(A/m)j cos A0J - Aj sin A0j 


(10) 


( 11 ) 
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Because this last equation is non-linear in the coefficients A< , it would cause 
difficulty in obtaining a solution for them. However, as mentioned previously 
only 5 of 6 coefficients in a given segment are independent; therefore, an addi- 
tional constraint can be imposed without affecting the resulting equation for 
the general conic section. The additional constraint used is to equate sepa- 
rately the numerators and denominators of both sides of Eq. (ll). This pro- 
cediire yields two equations, both linear in the coefficient Aj , which may be 
rearranged to give 


(A/n)j = 


jAj cos 


/sin A0 


J+1 


( 12 ) 


and 


(A/m)j = 


6-1 - 


/sin A0 


J 


(13) 


For the first segment (j = l),m2^ is specified and Eq. (13) is not needed, eind 
also for the last segment (j = N),nj^ is specified and Eq. (12) is not needed. 


Now substitute Eqs. (2), (3), (T)» (8), (l2), and (l3) into Eq. (l) to ob- 
tain the equation for the general conic sections as 


where 


“/j-i * ®/j * ''/j+i * 


= 0 


(IM 


a . = 




/sin A0 for J > 1 
J 

for J = 1 


(15) 


2 

y - (l/m^ + cotA02) yz + for J = 1 


3j = < y + (cotA0j - cotA0j+i) yz - cotA0jy^z - y^y for 1 < J < 


(16) 


2 

y + (cotA0jj - 1/njj) yz - cotA0jjy^jZ - y^^y for J = N 


Ya = 


yz/sinA0j^j^ for J < N 
for J = N 


(IT) 


These equations indicate that there are only two unknowns in each segment, Aj 
and Cj. However, there is one more segment than interior control points which 
means that one of these coefficients is strbitrary. Here Aq^ = 1 is used arbi- 
trarily (unless the conic section requires A^^ = O) and the independent unknowns 
become Cj , Aj(j = 2, •••, N). If there were as meiny data points as un- 
knowns, and if some conic section could be made to pass through these data 
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points, then Eq. (l^t) could he applied to all the data points to give (2N-1) 
linear equations for the (2N-1) coefficients C^, Cj , Aj for J = 2, •••, N. 
However, there are generally more than (2N-1) data points, and if Eq. (lU) were 
applied to all of them an overdetermined system of linear equations (See ref. 

8) would result. Therefore, a least-squares solution of the overdetermined 
system is used to detennine the coefficients C^, Cj, A< for J = 2, •••, N. 

This procedxire is described in Appendix B. After obtaining the coefficients in 
this manner, the inequality of Eq. (9) is checked for the possibility of com- 
plex roots for z in each segment. If the inequality of Eq. (9) is not satis- 
fied, the coefficient Cj in that section is replaced by the Value obtained us- 
ing the equality sign in Eq. (9). As mentioned earlier, the coefficient Cj 
does not affect the slopes of the ciirve at the end points of the segment, and 
hence C< does not affect other segments. The eqviality sign in Eq. (9) gives a 
hyperbola of 2 straight lines as shown in figure U. 

In order to solve for z as a function of y from Eq. (1I+), a quadratic 
equation must be solved, and the proper choice of the + or - sign must be de- 
termined beforehand for each segment. It is shown in Appendix A that in order 
to make z = 0 at y = 0 (a control point at the beginning of segment j), the 
+ sign must be used if (A/m)j > 0 and the - sign must be used if (A/m)j < 0. 

Once the coefficients Aj and Cj are determined, all the conic sections for 
that cross section are completely defined. In order to put these resvilts into 
a form suitable for "blending" the cross sections in the longitudinal direction, 
the conic sections for each segment are redefined in terms of 1* points: the 

two control points at the end of the segment, a slope point which determines 
the slope at the end points, and finally an intermediate point on the curve be- 
tween the end points (See figure 6). The 3 points on the curve and the two 
slopes at the end points of a segment are sufficient to determine new coeffi- 
cients A^, A2, A3, Ai|, and A5 for the general conic section*, 

—2 —— —2 — — 

A^y + A^yz + A^z + A^y + A^z +1=0 (18) 

in global coordinates y,z. This process is applied to each segment in a cross 
section, and therefore the 5 coefficients become functions of the longitudinal 
coordinate x when the segments of a cross section are blended with correspond- 
ing segments in the other cross sections. 


Longitudinal Variation of Cross Sections 

In order to determine the longitudinal variation of the coefficients in 
Eq. (18), a three-dimensional curve is fit through each of the 4 points used to 
define the conic section of corresponding segments (See figure 7). In contrast 
to the cross-sectional data points, these c\irves must pass through each of the 
points in the longitudinal direction. They are represented by their projections 


Note that the conic section given by Eq. (18) must have the constant 1 re- 
placed by 0 if the curve is to pass through the origin. 
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in the x, y and x, z planes; hence, two planar curves are used to represent 
each three-dimensional curve. The parameteric method of cubic splines (ref. 9) 
is used to curve-fit each planar curve, with the chordal distance between the 
coordinate points as the parameter. The parametric spline allows infinite 
slopes whereas the regular spline does not . 

This method worked very well for some test cases but had trouble with 
others because of the bumps and wiggles inherent to cubic splines. These bxomps 
and wiggles can be corrected by specifying slopes in the longitudinal direction 
at certain cross sections. Also parts of the longitudinal curves can be speci- 
fied as straight lines. 

Consider now the longitudinal variation of a conic section. In each cross- 
sectional plane, Eq. (l8) will hold but the coefficients A 2 , A 2 , A|^, and A^ 
will vaiy with x. As mentioned previously, these coefficients are determined 
by 4 defining points (the two control points, an intermediate point, and the 
slope point). For each segment, the 5 equations used to determine the coeffi- 
cients Aq (q = 1, •••, 5 ) are formed by applying Eq. (I 8 ) to the 3 points on 
the cross-sectional curve (the two control points z^ and jr^, Zj_, the inter- 
mediate point yjj, slopes at the ends of the segment using the slope 

point y^ , Zg (See figure 6). This procedure yields the following 5 equations: 


5 

1 


q=l 


G 

PI 


A 

<1 


d 

P 


P = 1. 


(19) 


At any longitudinal position, Eq. (l9) can be solved by any standard matrix in- 
version routine to determine the coefficients Aq. The derivatives dAq/dx and 
d^Aq^/dx^ can be obtained by differentiating Eq. (19) and successively solving 
the resulting system of linear equations. The elements of Gpq and their deri- 
vatives with respect to x are obtained from the three-dimensional curves which 
were spline-fit through the 4 points used to define the conic section for that 
segment in each cross-sectional plane. 


COMPUTATIONAL ALGORITHM 


Given the set of data points in cross-sectional planes at several 

longitudinal stations, 

( 1 ) For each cross-sectional plane, divide the data points into seg- 
ments so that a conic section can be curve-fit to the data points 
in each segment by the least-squares technique developed herein. 

( 2 ) If the curves fit to the cross-sectional data points are not 
satisfactory, modify them by one or more of the following methods: 
(a) define new boundaries (control points) for segments, (b) 
specify slope(s) at control point (s) (slopes may be finite or in- 
finite and continuous or discontinuous), (c) specify selected seg- 
ments as straight lines, (d) a specific conic section can be 
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specified for a segment by prescribing the slopes at the ends of 
the segment and using only one datum point between the end con- 
trol points . 

(3) Represent the conic section for each segment in a cross-sectional 
plane in terms of the two control points at the ends of the seg- 
ment, an intermediate point, and the slope point (See figure 6). 

{k) For each point foiind in step (3), spline- fit a three-dimensional 
longitudinal curve through it and the corresponding points in 
other cross-sectional planes (See figure T). 

(5) If the longitudinal curves are not satisfactory, modify them by 

one or more of the following methods: (a) specify slope(s) at 

longitudinal station(s) (slopes may be finite or infinite and 
continuous or discontinuous), (b) specify selected longitudinal 
segments as straight lines, (c) redefine the boundaries (control 
points) of the segments in the cross-sectional planes so that 
the points used for the longitudinal spline-fit form a smooth 
curve . 

(6) The geometrical properties of the surface may be computed in 
polar or Cartesian coordinates at any position x, y by the fol- 
lowing steps: 

a. Locate x between two consecutive longitudinal stations, 
and then locate the cross-sectional segment which con- 
tains y. 

b. Use the spline function to calculate the coordinates, 
slopes, and second derivatives of the 1+ longitudinal 
curves for this segment at x (See figure 7)* 

c. Calculate the coefficients A^(q = 1, •••, 5) of the 
conic section at this location by use of Eq. (19). 

Determine the first and second derivatives of Aq with 
respect to x from the first and second derivatives of 

Eq. (19). 

d. Calculate the body position z from Eq. (18), and the 
derivatives of z with respect to x and y from deriva- 
tives of Eq. (18) . 

(T) After a satisfactory surface fit has been obtained, the data which 
must be retained for a geometry subroutine package are the co- 
ordinates and longitudinal slopes of the longitudinal curves at 
those longitudinal stations where cross-sectional data points 
were given. Then, the geometrical properties of the surface can 
be calculated at any position by the method outlined in step (6) 
above . 
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APPLICATION TO 70° DELTA WING 


The siirface fitting method was applied to the 70° slab delta wing shown in 
figure 8. This example was chosen because it illustrates many of the options 
available to modify the longitudinal curves and because the results can be com- 
pared with an exact solution. 

Cross-sectional data was used as shown in figure 8, and due to symmetry 
only the first quandrant is used. Two segments (three control points) are 
needed to represent the cross section of x = 10. The first segment is a 
straight line and the second is one-fourth of an ellipse. The least-squares 
ciirve-fit technique represents the ellipse exactly by specifying a zero slope 
at control point j = 2, an infinite slope at control point j = 3, and one datum 
point between these two control points. Although two segments must also be 
used for the other two cross sections, only one segment is necessary to specify 
the circle at x = 0.65798 and the ellipse at x = 1.0. Therefore, the first two 
data points, which are also control points, are made coincident. Then the ex- 
act curves are calculated from the least-squares curve-fit by specifying a zero 
slope at control point j = 2, an infinite slope at control point j = 3, and one 
datum point between these two control points. The exact location of this datum 
point is irrelevant except that it must lie on the desired curve. 

The three-dimensional longitudinal curve is represented by its projections 
in the x,y and x,z planes. Therefore, l6 longitudinal planar curves are used 
for this example (8 for each cross-sectional segment). Modifications were made 
to the initial spline-fits to lU of these curves since exact slopes are readily 
obtained from figure 8. 

The geometrical properties and their derivatives were calculated at four 
circumferential positions for x = 0.3, 1.0, 2.0, and 5. The results at 
X = 1.0, 2.0, and 5.0 are exact (within the accuracy of single precision on the 
IBM 360/175 computer), whereas some inaccuracies were noted at x = 0.3. 


APPLICATION TO HL-IO LIFTING BODY 


The surface fitting method was also applied to part of the HL-10 lifting 
body. Figure 9 shows a view of this body and the six cross sections used to 
generate the surface fit. Tabulated values of the coordinates of siirface 
points were available for many cross sections in addition to the six actually 
used, and these additional cross sections were used as a check on the accuracy 
of the s\irface fit. Due to symmetry about the x,y plane, only half of the body 
was used. 

Figiire 10 illustrates the first attempt to curve fit the cross-sectional 
data at X = 128.32. Only two segments were used, and hence the resulting ciirve- 
fit is not satisfactory. Figure 11 shows the revised ciirve-fit to this same 
cross section when three segments are used, and the results are very good. 
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After satisfactory curve-fits were obtained for all six cross sections, 
using three segments in each one, the cross sections were blended in the longi- 
tudinal direction by use of the parametric cubic splines. Since each three- 
dimensional longitudinal cxirve is represented by its projections in the x,y and 
x,z planes, a total of 2k planar curves were used to represent the four longi- 
tudinal curves required for each of the three segments in a cross section. Of 
these 2k curves, seven of the initial spline-fits were found to have wiggles, 
and hence they were modified by specifying their slopes at the nose (X « O). 
Some of the revised longitudinal curves are shown in figure 9* and they are 
actually smooth althovigh the plotter makes them appear to have wiggles. 

As a check on the accuracy of the final surface fit, cross sections were 
calculated at axial stations in between those used to generate the surface fit. 
Figure 12 compares the calculated cross section at 5 • 43.6^6 with the actual 
coordinates for this body. Other cross sections were found to ccxnpare ' even 
more closely than this one with the actual coordinates. 


CONCLUDING REMARKS 


An algorithm has been developed for surface fitting three-dimensional bod- 
ies from data points in several cross-sectional planes. This method was found 
to give satisfactory surface fits to the 70® slab delta wing and HL-10 lifting 
body. The surface geometry can best be analyzed through the use of an inter- 
active computer graphics environment. 

1 

After a satisfactory surface fit to a body has been obtained, a relatively 
simple geometry subroutine package can be formed to use in other computer pro- 
grams requiring a mathematical model of the geometry. It will calculate the 
body coordinates, slopes, and second derivatives at any position on the body. 
This method could be vised to generate the input data for Coons' pathing method 
if desired. Computational time is short, and the amount of storage required is 
relatively small. A copy of the computer program can be obtained from the 
authors . 
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APPENDIX A 


Solution for Conic Section Equation 

Substitute Eqs. (2), (3), (7), and (8) into Eq. (l) to obtain the conic 
section in the local coordinate system in the form 


C z + P z + Q = 0 
J J J 


(Al) 


where 




(A/m) + (A/n), y + (A/m) y 


J 


J 


J •'J 


(A2) 


and 


Qj = Ajyjy-y 


The solution of Eq. (Al) is 


z = 


- P ± R 

J J 


/2C, 


^ 


where the discriminate (Rj) is given by 


'Jj 


for C 0 
J 

for C = 0, P yi 0 
J J 


(A3) 


(aU) 


■'j = '■j - ‘‘“jS • 


(A5) 


In order to obtain real roots for z from Eq. (aU), Eq. (A5) must give Rj ^ 0 
for 0 f y f yj* Note that Eq. (A5) gives Rj > 0 when AjCj > 0. 


Also, 


and 


Rj = (A/m)j^ >0 at y = 0 


Rj = (A/n)j^ >0 at y = yj . 


The possible variations for Rj are illustrated on the following page, 
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R 


J 


> 0, real roots 



The minimum value of AjCj vhich still gives real roots corresponds to the limit- 
ing case shown below. 



The minimum value of Rj occurs at y = y* where both R< = 0 and dRj/dy = 0. The 
solution of Ri = 0 and dR</dy = 0, using Eq. (A5), gives the minimum value of 


AjCj as 




= (A/m)j (A/n)j 


(a6) 


Thus, Eq. (A4) will give real roots for z if 


AjCj > (A/m)j (A/n)j . 


(AT) 


The choice of the sign to be used in Eq. (a 4) is determined by requiring the 
equation to be satisfied at the control points (z = 0, y = O) and (z = 0, 
y = yj ) • For the first control point, y = 0, Eq. (a 4) yields 


z 




2C 


(A/m)^y 

J 



(A8) 
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and th\is z = 0 requires the + sign if (A/m)j > 0 and the - sign if (A/m)j < 0. 
At the second control point, y = yj , Eq. (A4) yields 


z 





(A9) 


and z = 0 requires the - sign if (A/n)j > 0 and the + sign if (A/n)j < 0. 

These conditions are all compatible because a conic section passing through the 
two control points and single-valued in z will have nj < 0 if mj > 0 and nj > 0 
if mj < 0. The sign given by these conditions can be used in Eq. {kk) for all 
values of y in the range 0 < y < yj . 


APPENDIX B 


Least-Squares Solution for Conic Section Coefficients 

Equation (l4) cannot be satisfied at all data points because there would 
be more equations than unknowns. It is also not desirable, in general, to have 
the curve go through all the data points because there may be scatter in the 
data. Accordingly, Eq. (14) is applied at data point k and rewritten as 


“j,k ^J-1 


.k 


^J,k ^j+1 




■J.k 


(Bl) 


where Yj k "^ke values of oj , 3j, and yj evaluated at 

y = yj^, z = Zk*i^ segmen"^ J; and rj ^ is called the residual. The least-squares 
solution of the overdetermined system of equations determines the coefficients 
Cj, (^ , Aj (j = 2, •••, N) which minimize the sum of the residuals squared (See 
ref. 8) . Define Kj as the data point n\imber which corresponds to the first 
control point in segment J. Square Eq. (Bl) and sum over all the data points 
in continuous segments J = 1, •••, N to obtain 



The right side of Eq. (B2) is minimized by the system of equations obtained by 
setting partial derivatives of Eq. (B2) with respect to the independent coeffi- 
cients equal to zero. The resxxlt of setting partial derivatives with respect 
to Ap, • • • , Ajj equal to zero yields the following set of equations: 
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The result of setting partial derivatives with respect to C;j_, •••, equal to 
zero gives the following additional system of equations: 


h - + b -A + b oA I C =0 

m,l m-1 m,2 m m,3 m+1 m,4 m 

m = 1, • • • , N . 


(bM 


The combined system of Eqs. (B3) and (bU) gives a system of (2N - 1) linear 
equations for (2N - l) coefficients. The peurameters used in Eqs. (B3) and (bU) 
are defined as follows: 
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Note that if slopes should 

be specified 

at selected control points. 

then 


the least-squares solution described above is applied to the segments between 
two consecutive control points with specified slopes. In the solution of the 
combined system of Eqs. (B3) and (B4) the term Ufj+j and Yfj+i^k inust be inter- 
preted as zero. Hence, the coefficients A^, and Aj ^^.2 do not appear in 
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the resiilting system, and recall that = 1 unless the conic section requires 

= 0 . 
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Figure 1. Control points and data points in a cross section. 
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Eqiiatlon for General Conic Section 


Aj + Bj yz + Cj y + z + Fj * 0 


Z 
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O CONTROL POINT, and (y^ ) 

O INTERMEDIATE POINT, (y^.i^) 

□ SLOPE POINT, (yj.ij) 



Figure 7. Longitudinal curves through the four points used to define 
a segment of a conic section. 
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x=10 (straight line + ellipse) 


Figure 8. Geometry of 70° delta wing. 
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Figure 9. View of cross sections and longitudinal curves for HL-10 body. 
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